Perturbation Method for Classical Spinning Particle Motion: II. Vaidya Space-Time 
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This paper describes an application of the Mathisson-Papapetrou-Dixon (MPD) equations in 
analytic perturbation form to the case of circular motion around a radially accreting or radiating 
black hole described by the Vaidya metric. Based on the formalism presented earlier, this paper 
explores the effects of mass accretion or loss of the central body on the overall dynamics of the 
orbiting spinning particle. This includes changes to its squared mass and spin magnitude due 
to the classical analog of radiative corrections from spin-curvature coupling. Various quantitative 
consequences are explored when considering orbital motion near the black hole's event horizon. An 
analysis on the orbital stability properties due to spin-curvature interactions is examined briefly, 
with conclusions in general agreement with previous work performed for the case of circular motion 
around a Kerr black hole. 

PACS numbers: 04.20.Cv, 04.25.-g, 04.70.Bw 



I. INTRODUCTION 

The Mathisson-Papapetrou-Dixon (MPD) equations 
[1-4] represent a well-known description of classical spin- 
ning particle motion in the presence of a curved space- 
time background. They comprise the "pole-dipole ap- 
proximation" for the dynamics of extended bodies with 
spin angular momentum in the vicinity of black holes, 
neutron stars, or other sources of space-time curvature 
where a strong gravitational field is generated. This in- 
cludes sources which themselves are time- varying for the 
duration of a spinning particle's motion along its world- 
line. As a consequence, the spin-curvature coupling term 
in the MPD equations, which generates an external force 
and torque to act on the spinning particle, also becomes 
time-varying, leading to potentially very interesting dy- 
namical effects experienced by the particle. 

One particularly interesting space-time background 
with an explicit time dependence is known as the Vaidya 
metric, which describes space-time curvature due to a 
spherically symmetric compact source that either radi- 
ally accretes surrounding radiation, or radiates away its 
central mass. Given that most astrophysical sources have 
at least some orbital or spin angular momentum during 
their formation, it is unlikely to find candidate sources 
in the night sky that carry the properties exactly de- 
scribed by the Vaidya metric. However, because of its 
relative simplicity compared to the Kerr metric to de- 
scribe rotating black holes, while also having a time- 
dependent central mass, the Vaidya metric nonetheless 
provides an ideal testing ground for understanding sub- 
tle properties of the MPD equations for an orbiting spin- 
ning particle that is sensitive to a time-varying gravita- 
tional field. A recent paper [5] presents an extensive nu- 
merical investigation of the MPD equations in a Vaidya 
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background, modelling a point dipole in circular orbit 
around a much heavier non-rotating black hole described 
by a monotonically increasing central mass function in 
terms of known functions. This paper also shows, using 
only the quadrupole moment formula, that the dynami- 
cal background due to a growing central mass can influ- 
ence the shape and frequency of gravitational waveforms 
generated by the spinning particle for a sufficiently large 
mass accretion rate, with potentially useful implications 
for low-frequency gravitational wave astronomy via the 
space-based LISA observatory [6]. 

Although a numerical treatment of the MPD equations 
in a Vaidya background is undoubtedly a useful exercise, 
an analytical exploration of the same problem is defi- 
nitely beneficial in many respects. For example, knowing 
the explicit time-dependence of the mass function within 
an analytical expression of the MPD equations allows for 
the study of conditions where instabilities in the dynam- 
ical system most likely will occur. It can also potentially 
give useful insight for knowing when a mass increase or 
loss will lead to macroscopic changes in the particle's or- 
bit for a predetermined mass accretion or loss rate. Fur- 
thermore, the general results obtained from such a study 
can provide clues for how a spinning particle may respond 
due to a more realistic time-dependent source than one 
described by the Vaidya metric, such as a pulsating star, 
particularly on determining the most dominant contribu- 
tion to its response. 

A recent development on the study of the MPD equa- 
tions involves a linear perturbative approach first intro- 
duced by Chicone, Mashhoon, and Punsly (CMP) [7], 
with an application by Mashhoon and Singh [8] for de- 
termining a first-order perturbation of a circular orbit 
around a Kerr black hole due to spin-curvature cou- 
pling. This first approach was more recently generalized 
by Singh [9] to accommodate for higher-order contribu- 
tions in powers of s/(mr) <^ 1, where p = s/m is the 
M0ller radius [8, 10] in terms of the particle's spin mag- 
nitude s and mass m, and r is the radial distance from the 
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background mass source to the particle's location. This 
generalization can be applied to formally infinite order in 
the perturbation expansion parameter and makes no ref- 
erence to any particular space-time metric or symmetries 
therein. A detailed application of the generalized CMP 
approximation to the MPD equations was just presented 
for the case of circular motion around a Kerr black hole. 
For future reference, this recent paper is now identified 
as "Paper I" [11]. It would be very interesting to perform 
the same investigation as found in Paper I, but this time 
applied to the Vaidya metric. 

The purpose of this paper is to apply the generalized 
CMP approximation of the MPD equations to describe 
circular motion around a static compact object described 
by the Vaidya metric, and incorporate both mass accre- 
tion from null radiation and outgoing radiation within 
the formalism. This paper begins with a brief review of 
the MPD equations and the generalized CMP approxi- 
mation [8, 9, 11], found in Sec. II. An introduction to 
the Vaidya metric [5] and its application to the gener- 
alized CMP approximation is then presented in Sec. III. 
Following this. Sec. IV describes the main results for the 
case of circular motion around the central body to sec- 
ond order in the perturbation expansion parameter, in- 
cluding the "radiative corrections" of the squared mass 
and spin magnitudes predicted within the underlying for- 
malism [9, 11]. Afterwards, a discussion of the obtained 
results is given in Sec. V, followed by a brief conclusion. 
Consistent with Paper I, the Riemann and Ricci tensors 
follow the conventions of MTW [12] with signature +2, 
and assuming geometric units of G = c = 1. 



The MPD equations presented in (1) are underdeter- 
mined, and require supplementary equations to specify 
the system. A commonly accepted constraint is to im- 
pose orthogonality between the particle's linear and spin 
angular momenta, following Dixon's approach [3, 4], such 
that 



S^'^Pg = 0. 



(2) 



As well, the mass and spin parameters m and s are iden- 
tified by the constraint equations 



m 



(3a) 



(3b) 



which become constants of the motion [7] when (2) is 
implemented within the MPD equations. It is also well- 
known that the four- velocity can be expressed in terms 
of Pi^ and 5"^ within the MPD formahsm [13], leading 
to 



P- 



1 g"^-^ R^^^o^p P-' 



4 Raflpa S"'^ S^"^ 



, (4) 



where specification oi P ■ u determines the parametriza- 
tion constraint for r. Clearly, (4) shows that the spin- 
curvature coupling creates a displacement of the parti- 
cle's four-velocity away from geodesic motion. 



B. Generalized CMP Approximation 



II. MATHISSON-PAPAPETROU-DIXON (MPD) 
EQUATIONS AND THE GENERALIZED CMP 
APPROXIMATION 

A. MPD Equations 

Given the dynamical degrees of freedom P''(t) and 

S"^{t) for the spinning particle's linear four- momentum 
and spin tensor, respectively, the MPD equations are 



DPt^ 
dr 

dr 



(la) 
(lb) 



where R^^vafi is the 
u''(t) = da;^(r)/dT 
parametrization r. While 



Riemann curvature tensor and 
is the four-velocity with affine 



T can be chosen to satisfy 



w'^^i/i = —1 to describe proper time, it is not necessary 
to impose this particular constraint if desired. The com- 
bined force equation (la) and torque equation (lb) infer 
that the particle's four-momentum precesses around the 
centre-of-mass worldline, giving rise to non-trivial motion 
away from time-like geodesic motion. 



While a more detailed account of the generalized CMP 
approximation can be found in Paper I, it is useful to 
briefly summarize the main points of this approach to 
the MPD equations. This is a perturbation approach 
based on the assumption that 



(5a) 



3=0 



S'^^e) ^ sY^e^S^^ = J^e^S^;_^^, (5b) 



where P^^^ and S'^j^i^ are the respective jth-order con- 
tributions of the linear momentum and spin angular mo- 
mentum in e, an expansion parameter associated with s. 
In addition, the four-velocity is described as 



\3) ■ 



(6) 



The zcroth-order expressions in e then correspond to a 
spinless particle in geodesic motion, while higher-order 
contributions are identified with spin-curvature coupling. 

The main idea to the generalized CMP approximation 
is to substitute (5) and (6) into both the MPD equations 
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(1) and the exact expression for according to (4), ex- 
pand these equations with respect to e, and solve for each 
order of the perturbation expansion iteratively. It follows 
that the jth-order expressions of the MPD equations are 



W 

dr 



0-1) 
dr 



2 



fe=0 



fe=0 

where j = implies that 
dr 



-k) "(fe) ' 



0, 



(7b) 



(8) 



while j — 1, corresponding to the CMP approximation 
[7, 9, 11], is 



DPm 1 R 

~ 2 ''"^ (0) (0) ' 



DS: 



(0) 



dr 



= 0. 



(9a) 



(9b) 



In addition to (1), the supplementary spin condition 
equation (2) and constraint equations (3) for the squared 
mass and spin magnitudes need to be incorporated within 
this formalism. For the spin condition, it is straightfor- 
ward to show that 

^r^[;) = -E^f c-fc)' i>i (10) 



fc=i 



for the (j-l-l)th-order contribution, where 



(11) 



for the first-order perturbation in £. As for the squared 

mass and spin magnitude constraint equations, it is possi- 
ble to identify a bare mass nio and bare spin sq according 



to 



si 



_p(0) 

1 0(0) Qtil' 

2 '^tif "^(0) ' 



(12a) 
(12b) 



such that 



m^(ff) = mg I 1 + E^ 



where 



1 Ypl^r'^pr 



m. 



2 ^ n 

k=0 



(k) 



1 ^ 

„2 '^tii^ ^ik) ' 



(13a) 
(13b) 

(14a) 
(14b) 



fc=0 



are dimcnsionless jth-order "radiative corrections" to 
TOq and Sg, respectively, due to spin-curvature coupling. 
Each expression of (14) satisfies 



dr 



Drrii 



(15) 



Solving for the four-velocity (6) requires specifying the 
parametrization constraint within (4). Following Paper I, 
the particularly useful choice of 



P-u 



-m(e). 



(16) 



leads to 



4 



3=0 



U) 



mo 



+ £ 



mo 



^(1) 



+ £M — 
Too 



Pi 



(2) 2 ^ (1) 



3 -4 



^(0) 



a/3 



2to3 (0) -f^i'7a/3 -f^(o) -^(O) 



Too 



-^(3) 2 ^ (2) 



to! - ^m^ ) Pj'^) - ^ (mi -:Linlml+^-inl\P(^^ 



+ 



2ml 



EqUiy p7 r,a;/3 ™2 qMI^ p7 C 



Ln=0 



fc=0 



(17) 



where 



(18) 



implying that li'' is indeed the four-velocity with unit 
normal to third-order in e. It is not necessarily true, 
however, that u-u = —1 applied to all orders of e corre- 
sponds to (16). Extending (18) to fourth-order in e and 
higher requires a more general approach, where 



P -u 



3=0 



(19) 



and the choice for each (P • u)^j^ is determined from con- 
straint equations for each order of s as required. 



C. Summary of the Linear Momentum and Spin 
Angular Momentum Expansion Components 

The approach adopted to solve for the linear momen- 
tum and spin tensor expansion components in the gener- 
alized CMP approximation is to use the tetrad formalism 
and work in Fermi normal co-ordinates. Full details for 
obtaining these expressions are shown in Paper I, but 
it is worthwhile to give a brief outline of the procedure. 
Suppose that an orthonormal tetrad frame A'*^ satisfying 



V. 



&0 



and parallel transport (-DA'^a/dr = 0) describes a projec- 
tion of space-time curvature described by general space- 
time co-ordinates /i onto a locally flat tangent space, de- 
noted by a. The Fermi co-ordinates are described by 
X" in the local neighbourhood about the spinning par- 
ticle's centre-of-mass worldline, while general space-time 
co-ordinates are denoted by Jf*. As usual, A^g = u^^^y 
Furthermore, the Riemann curvature tensor in the Fermi 
frame is then given by 



■R 



afj'yfi ~ ^IJ-vpiy ^'^a A X'' j A j , 



(21) 



and that for j >0, 



(22a) 
(22b) 



Following the approach taken in Paper I and elsewhere 
[9, 11], it is shown from (22) for j = that 



pM 

Uo) 

•5(0) 



A'^Q P(o) — 



mo A'*e , 



(23a) 
(23b) 



satisfying the first-order spin condition (11), where 
P^p^ — Too (5"q and S'^^^ is a constant-valued spatial anti- 
symmetric tensor determined from initial conditions. 

For J = 1, the linear momentum is straightforwardly 
determined to be 



while 



-^(1) 



^(1) - 2 



\IJ. \1/ Q&fS 



(24) 



4 ^ mo J I 



'(0) 



(20) subject to 



dr 



0. 



.(25) 



(26) 



Contracting (24) into P^°^ shows that the first-order mass 
shift contribution is identically 



TOi 



0. 



(27) 



However, the expression for first-order spin shift is still 

formally undetermined based on (25) alone, and while it 
is tempting to set sf = in analogy with (27), this is 
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not justified given that only needs to be covariantly 
constant according to (15), and not necessarily zero. To 
obtain an expression for requires the direct solving of 
(26), the details of which are given in Paper I and are 



presented in Appendix A of this paper. 

Solving for the j = 2 expressions for both and S"^ 
is straightforward, such that 



for the linear momentum, where 



= 1^ a1(r)dr 



(29) 



is the time-averaged jth-order correction to the squared 
spin magnitude. The corresponding expression for the 
spin tensor is given by 



(30) 



the solution to 



^^i"^) _ J_ p[M e""]"" E> . p7 qa(3 /oi\ 
dr ~ ^3 -^(0) '^(0) -^<^7a/3 -^(0) '^(0) ' V"*"^/ 

after substituting u'^^^ from (17). 



D. Perturbations of the M0ller Radius 

As noted in Paper I, the M0ller radius p = s/m is 
closely identified with the strength of spin-curvature cou- 
pling experienced by the spinning particle. Previous 
studies of chaotic dynamics in the Kerr background [14- 
17] indicate the possibility that perturbations of p may 
reveal the conditions where a transition from stable to 
chaotic motion can appear for a spinning particle in a 
general space-time background. The perturbation ex- 
pression for the M0ller radius is then formally given by 



m{e) 



Too 



1 



(^1 



to: 



+ 0(£3)}, 



(32) 



where the contributions due to rh\ are retained for com- 
pleteness' sake. It is of particular interest to see how 

(32) behaves for the Vaidya metric, which can then be 
compared directly with the results obtained in Paper I. 



III. GENERALIZED CMP APPROXIMATION 
IN VAIDYA SPACE-TIME 

With the formalism of the generalized CMP approxi- 
mation presented, it is possible to now develop the frame- 
work for applications to motion in a Vaidya space-time 
background. The most immediate challenge is to derive 
the orthonormal tetrad frame A'^a for application of the 
formalism just outlined. It is very surprising to note that, 
while the Vaidya metric is much simpler in form com- 
pared to the Kerr metric used in Paper I, the relevant 
computations are technically much more involved, lead- 
ing to much greater complexity than first anticipated. 
This is because the Vaidya metric is effectively time- 
dependent, since the mass function is no longer static, but 
either grows or shrinks monotonically along null rays. Ul- 
timately, this property must be incorporated within the 
structure of the orthonormal tetrad. 

The Vaidya metric in {^,r,6,(j)) co-ordinates is de- 
scribed in general form as [5, 18, 19] 



ds^ 



1 



d^^ -|-2ad^dr 



+ (d6l^ -F sin^ Q ( 



(33) 



where ^ is a generalized null co-ordinate denoting time 
development and a is a dimensionless parameter chosen 
such that 



a 



(34) 



for ingoing radiation [5] along the advanced null co- 
ordinate while 



a. = —1 



(35) 



corresponding to outgoing radiation [19] along the re- 
tarded null co-ordinate /it. For the Vaidya metric, the 
central mass function Af(^) is a monotonically increas- 
ing or decreasing function of ^(r) for a given choice of 
a to satisfy the weak energy condition, but is otherwise 
an arbitrary function. The mass function can also be 
defined as 



M(0 = Mo + AM(e) 



(36) 
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where Mq is the static mass for a Schwarzschild black 
hole and AM(0) 0. It will prove useful to express (33) 
in terms of {t, r, 6, (p) co-ordinates, where ^ is described 
by the tortoise co-ordinate condition [5] 



t + a 



r + 2Mo In 



2Mn 



1 



(37) 



This leads to the Vaidya metric expressed as 



ds^ = - 



2Mr 



1 - 



4a 1 - 



2 AM 



2Mn 



dtdr 



+ 



AM 



+ r2 (d^2^sin2^d<^2^ , 



(38) 



which reduces to the Schwarzschild metric as 
AM(^) ^ 0. 

While it is mathematically acceptable to leave AM(^) 
unspecified, it creates computational obstacles for an ex- 
act treatment of the problem. Therefore, a simplifying 
assumption adopted is to let AM/Mq <C 1, which is well- 
justified on physical grounds, since the Eddington lumi- 
nosity limit [20] imposes an upper bound mass accretion 
rate of 

where Mq is one solar mass and 7 « 0.1 is the energy 
release efficiency of the outgoing photon flux. This allows 
for a derivation of the Vaidya tetrad frame A^^ in terms 
of a linear perturbation about A^q, (Sch)) the Schwarzchild 
tetrad frame for circular motion, which is presented be- 
low. 

Consider the Schwarzschild orthonormal tetrad frame 
[8] for circular motion with fixed radius r > 2Mo, such 
that 



'^''o (Sch) 

'^''i (Sch) 



E 



,0,0 



L 



*(Sch)' 



A^' ' 'r2 sin6i 

— ^ sin (r^K t) , A cos (fix t) , 0, 
r A 



(40a) 



E 



tA sm6 



sin (Qk ' 



A'' 



2 (Sch) 



1 



0,0, -,0 
r 



A" 



3 (Sch) 



— 7 cos (Qk t) , A sin (Qk t) , 0, 

rA 



E 



r A smO 



cos (f^K r) 



where 



Ok = 



(40b) 

(40c) 

(40d) 
(41) 



is the Keplerian frequency of the orbit, 



N 



3Mn 



2Mn 



(42) 



and the energy E and orbital angular momentum L for 
the orbit are 



E = — 



4! 

N 



rHh 



(43) 



The boundary conditions are determined such that t = 
= at r = 0. Furthermore, given that the Vaidya met- 
ric is spherically symmetric, the Cartesian axis centred 
on the black hole is oriented such that orbital motion is 
confined to the plane defined by ^ = 7r/2 with respect to 
an assigned z-axis. 

The next step is to derive the Vaidya orthonormal 
tetrad in the form 



A". 



A'*, 



(Sch) 



+ AA''^ 



(44) 



where AA'^a is the linear perturbation proportional to 
AM. While an exact treatment within this perturbation 
approach is given in Appendices B and C, the outcome 
is considerably more complicated than for the exact or- 
thonormal tetrad in the Kerr background [8] . Therefore, 
another simplifying assumption is introduced, in the form 
of a series expansion with respect to inverse powers of N, 
since any deviations away from the Schwarzschild con- 
tribution will only be potentially identifiable when the 
spinning particle approaches the nearest (photon) orbit 
of r — > 3A'/o, corresponding to ^ 0. In addition, an 
expression for AM needs to be chosen in terms of r that 
is consistent with both the properties of the metric (38) 
and the mass accretion rate upper bound (39). This leads 
to the choice of 



AM(t) 



a 
A 



d(AM) 



(45) 



where the prefactor of a in (45) accounts for the direction 
of radiation fiow, and |d(AM)/d^| < 1. With these 
further assumptions incorporated, it can be shown that 
the Vaidya orthonormal tetrad frame components for a 
particle in orbit near the event horizon {N 0) a,t 8 = 
7r/2 are 



1 



A^ « ^ 



N 2iV7 



d(AM) 



(rnK)C(r,OKr) , 

(46a) 







2a 



d(AM) 



sin^ (fix t) 



\2 











0, 



(46b) 

(46c) 
(46d) 
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for A'' 



0' 



1 



_NA 2Af7 
X sin (JIk , 



d(AM) 



2a 
7V3 



d(AM) 



A^i = 0, 



for A", 



X sin (Qk t) , 



A°2 = 0, 
A^- - ^ 



d(AM) 



d^ 



C(r, 17k r) 

(47a) 

sin^ {Qk t) , (47b) 
(47c) 
{rnK)Cir,ftKT) 

(47d) 



for A'' 2, and 



A°6 » rOK 



1 a 
+ 



X\ = 0, 
A's = 0, 



d(AM) 



(48a) 
(48b) 



NA 2m 



d^ 



C(r,OKT) 



X COS (r^K r) 
A^g » ^ sin(f7KT) 



(49a) 



2a 
7V3 



d(AM) 



\\ = 0, 



dC 

X COS (f^K t) , 



sin^ (Ok t) cos (Ok r) , (49b) 



d(AM) 



d€ 



(49c) 

(49d) 



for A'^g, where 



C (r, Ok t) 



2 sin (2 Okt) 

[(1 - 2 r Ok) sin (2 Ok r) - 2 Ok r] . 



r Ok 



(50) 



Given (46) (49), it is now possible to obtain the Rie- 
mann tensor components in the Fermi frame. While the 
exact expressions for ^-R^£>g^ are found in Appendix D, 
for the special case of ^ and 6 = tt/2 considered in 
this paper, the dominant nonzero components are 



■0101 



■0103 



■0113 



0212 



02 

"Ar2 

F 



2A^ + r^ O^ + 
2 ^2 + 0| + 



3 a 
3 a 



d(AM) 



dC 

d(AM) 



d^ 



(r^ O^) C(r, Ok r) 
(r^ 03,) C(r, Ok r) 



cos2(Okt) = --^iiasSs' 



sin (Ok t) cos (Ok t) 



223 ' 



3 



K 



iV2 

o2 



A(rOK) + 



a 



-"-0202 ~ ^2 



3 a 



N6 

d(AM) 



d(AM) 



(r'^ 0^)C(r,OKr) 



dC 



30| 
" 7V2 



^(-^k) + ^ 



dC 

(r3 0|) C(r, Ok t) 
d(AM) 



cos (Ok t) = — ''^i?! 



0223 ' 



— 

~ -"-131 



313 ' 



dC 



(r^ 0^)C(r,OKr) 



sin (Ok t) 



0313 ' 



(51a) 

(51b) 
(51c) 

(51d) 

(51e) 



0303 



O^ 

___K 

-^R 



2 A' 



3 a 



d(AM) 
d€ 



(r^ 0|) C(r, Ok r) 



sin^ (Ok t) 



1212 • 



(51f) 



It is straightforward to confirm that (51) agrees with (55) of Paper I in the Schwarzschild limit when a = 0. As well, 
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it is interesting to note that the Vaidya contribution to where 
the curvature becomes significant when |d (AM) /d^| 
A^^, which sets an appropriate scale for the mass accre- 
tion or loss rate in the analysis to follow in this paper. 



IV. APPLICATION TO CIRCULAR MOTION IN 
THE VAIDYA BACKGROUND 

Having now obtained the orthonormal tetrad for the 
Vaidya background in the Fermi frame, it is possible to 
make use of the generalized CMP approximation for the 
MPD equations. This first requires evaluation of the un- 
perturbed orbit from P^Q^(r), which then gets integrated 
with respect to r to eventually obtain (r) . Because 
of the time-dependence in the Vaidya metric due to the 
evolving central mass function, it is clear that the unper- 
turbed orbit will not be truly circular. However, since 
|d(AM)/d^| <C 1, the deviation from circular motion 
is minimal. Based on (23a) and (46), the unperturbed 
four-momentum components are determined to be 



^(°o)(r) = 



4)M = 



mo 



1 a 



N 2m 



d(AM) 



r Ok C (r, Ok t) , 
(52a) 



a 



d (AAf) 



0, 

mo Ok 



sin^ (Ok t) , (52b) 
(52c) 



1 

TV 



+ 



a 



2 7V7 



d(AM) 



d^ 



rOKC(r, Okt) 



(52d) 



It is clear from (52b) that the radial component of the 
four-momentum is directed inwards for infalling radiation 

{a — V) and outwards for outflowing radiation [a — —1). 
This makes physical sense because a growing central mass 
creates stronger curvature that gives rise to a stronger 
inward force felt by the orbiting particle, and vice versa 
for a dissipating central mass. Integrating (52) over t in 
the form 

^ro) (-) = ^J^^ % (^') + ^io) (0) (53) 

leads to the unperturbed orbit, where the initial position 
is X^^(O) = (0, r, 7r/2, 0) to correspond with the ar-axis. 
Evaluation of (53) results in 



a 



N 2N'' 



d(AM) 



xr^K [ C{r, Ok r') dr' , (54a) 
Jo 



Ok iV3 



d(AM) 



dC 



X 

2 ' 



Ok t — - sin (2 Ok t) 



(54b) 
(54c) 



Xfo)ir) 



Ok T a Id (AM) 



N 2iV7 I d^ 

xrO|/ C (r, Ok r') dr' , (54d) 
.In 



f C(r,OKT')dT' = {2 r Ok [1- cos (2 Okt)] 

Jo Zrllj^ 



+ Ar[(l-2rOK)[l-cos(20KT)]-20^r2]} 



(55) 



It is clear from (54) and (55) that a spinlcss particle in 
the Vaidya background experiences a quasi-circular orbit 
with an overall growth or decay of its radial position over 
proper time, plus some non-trivial oscillatory structure 
embedded within its time development. 



First-Order Perturbations in e 



As with the computation in Paper I, progressing to 
the first-order perturbation (CMP approximation) in 

the Vaidya background is conceptually straightforward. 
However, the outcome is analytically more complicated 
than for its counterpart in the Kerr background. For the 
spinning particle initially positioned on the .T-axis of the 
Cartesian frame, the initial spin orientation (0, (j)) for S^^^ 

is chosen [8, 11] to agree with the standard definition of 
{9, (j>) for the spherical co-ordinates with respect to the 
Cartesian frame's z-axis. This leads to 



c23 
*(0) 

-^(0) 

ci2 
■^(0) 



with the outcome that 



So sin^ cos( 
— So cos^. 
So sin^ sin^ 



(56a) 
(56b) 
(56c) 
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(57b) 
(57c) 
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sin(2 Ok T + 61) + sin(2 Ok r - 61) - N cos9 (0; 



iKT)| 



(57e) 



(57f) 



As noted earlier in Paper I, a complicated beat structure 
in the sinusoidal functions exists in (57), due to the ini- 
tial spin orientation angles. It is also confirmed that the 
leading-order spin tensor agrees with its Paper I counter- 



part in the Schwarzschild limit when a = 0. 

The first-order perturbation of the linear momentum 
is determined from (24). A straightforward evaluation 
leads to the expressions 
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It is useful to note the ratio between the azimuthal com- 
ponent of the linear momentum to its time component, 



P(°i)(r) 
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T 



1 - 
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d(AM) 



d? 



X < sin(2 Ok r) — 



TV 
(t-Ok) 



(Okt) 



• (59) 



Comparison between (59) and its counterpart (61) in Pa- 
per I shows an important distinction between the two 
expressions, where; (59) records a predominantly oscilla- 
tory f,ime variation in the ratio, while the expression in 
Paper 1 gives a strictly time-independent ratio of E/L. 
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B. Higher-Order Perturbations in e 

While computations for the second-order perturbation 

quantities arc straightforward to perform, in the Vaidya 
background they become prohibitively long. Therefore, 
expressions for the second-order linear momentum and 
spin tensor components arc presented in numerical form, 
contained in the next section. In similar fashion to that 
in Paper I, the "radiative corrections" to the squared 
mass and spin magnitudes (14a) and (14b) exist in rela- 
tively compact form, necessary to evaluate the perturbed 
M0llcr radius (32). 

It is first necessary to introduce the notation for beat 
functions in the form 

{ill T , n2 , ns (j)) = cos{niD,K t + n2 - (j)) 

± cos(nir2K T — 71.2 6* — ns 0) , 

(60a) 

Qf{niflKT,n29,ns$) = sm{niQ.KT + n2 -n3 4>) 



± sm{niClK T — n2 — Tis t 



Then the first-order spin shift in the squared spin mag- 
nitude in the Vaidya background is 



So 
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(60b) where 
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Q+{0,30,2^)-Q+{0,§,2^) +2 cos(3 ^) - cos ^ | , 



(62a) 



^ (r* ni) [24 (f2| r^) {2 [g+ (4 r , 3 ^ , 2 .^) - Q+ (4 Ok T , ^ , 2 4>) 
+ 3 [Q+(2OKT,3^,0)-Q+(2OKr,^,0)] - Q+(0 , 3^ , 2 0) + Q+(0 , ^ , 2 0)} 

+ (0Kr){l2 [Q+(40Kr,3^,2(^) -Q+(40Kr,e^2, 

+ ^ (5A + 8rOK) [Q+{2nKT ,39 ,0) - Q+{2nKT ,e ,0) 

48 r , - , . 

+ — {A + GrflK) \Qt{0,i9,2(j))-Q+{0,9,2^ 

-J (15^-8rOK) [0+(4nKr,3^,2<^)-Q+(4OKT,^,2, 
-J {27A-328rClK) Qt{2flKT ,39 ,0) + j (9A-h232rOK) Qt{2ClKT ,9 ,0) 
+ 1 (SA-AOrftK) [Qti^KT ,39 ,24,) - Qti^KT ,39 ,0) - Qti^KT ,9 ,2^) + Q+i^KT ,9 ,0) 
+ ^ (3^ + 8rOK) [o+(0,3^,2<^)-g+(0,^,2<^)] . (62b) 



When integrated over a cycle defined by the Keplerian frequency, the time-averaged expression for (61) is deter- 
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mined with respect to 
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27r/nK 



^ (r^ O^) sin^ ^ cos^ 3 cos(2 ^) - 1 
sin^^'cos^ 

16 A 



(63a) 



It is clear from (63) that (sf ) is well-behaved for the full 
range of 9 and (p. 



- (21 A - 32 TT^ A + 32 r I^k) sin(2 ^) 

For the second-order mass shift in the squared mass 
12 7r (A + 12rf2K) cos(2 (^) + 367r^ .(63b) magnitude, it is shown that 
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(64) 



where the corresponding time-averaged expression is 

2 
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(65) 



In similar fashion, the second-order spin shift is deter- 
mined to be 
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whose time-averaged expression is 



3^2 / So 



(r' Q-^) sin^ cos 9 (2 cos^ (j) - 1) 



+ 



1 



So 
mor 



d (AAf) 



|84a^(r^O|:) sin2^(4cos2^-l) sin(^cos^+^ (sLsi6)| • (67) 



V. NUMERICAL ANALYSIS 

Based on the analytic expressions presented in this pa- 
per, it is uscfiil to explore some numerical analysis for the 
main results obtained. For the same reasons as given in 
Paper I, the purpose for taking this route is to visually 
determine the consequences of increasing the order of the 
perturbation expansion in the generalized CMP approxi- 
mation, when applied to the Vaidya background. It may 
be possible to identify a correspondence between this ap- 
proach and that of a purely numerical treatment of the 
MFD equations. This numerical treatment of the gen- 
eralized CMP approximation is given in terms of plots 
found in Appendix E. 

As with Paper I, the emphasis for this analysis is to 
identify the stability properties of the spinning particle's 
motion in the Vaidya background, assuming r — 6M, 
= (j) ^ 7r/4, and |d(AM)/dC| = IQ-*. While it is 
understood from (39) that this choice for the mass ac- 
cretion rate (a = 1) is too large given the Eddington 
luminosity limit, it nonetheless provides a useful means 
to directly compare with the corresponding set of plots 
found in Paper I. As for the mass loss rate (a = — 1), 
there is apparently no reason to exclude this choice for 
|d (AM) /d^l on physical grounds. 

Throughout this paper, the numerical analysis assumes 
that /Li = so/{mor) = 10~^ and /x = 10"""^, where mo = 
10~'^M. To determine the magnitude for a realistic spin, 
it is first shown that 

so = (lO^ ml (68) 
for given r and mo, which suggests that 

M<10-^V (69) 

to accommodate for a realistic spin of Sq ^ m^, cor- 
responding to solar mass black holes and neutron stars 
[16, 21] in orbit around supermassive black holes. Since 
/X = 10~2 and /x = 10~^ lead to unrealistically large val- 
ues [16, 17] for So, while agreeing with the choice given 
previously [14, 15] to explore chaotic behaviour for the 
MPD equations, it follows that any chaotic effects de- 
termined in this paper occur outside of astrophysically 
realistic conditions, as noted in Paper I. 



Figure 1 lists the plots of the particle's co-ordinate 
speed 

v{s) = ^gijV\e)Vi{e) (70) 
as a function of r, where 



for u'^{e) according to (17), to first- and second-order in e, 
and with the restriction of < w < 1 . As with the result 
of Paper I, Figs. 1(a) and 1(b) show that the expression 
to first-order in e is almost exclusively responsible for v 
due to So/ {mo r) = 10"^. Unlike the corresponding set of 
plots due to the Kerr background, however, the range of 
co-ordinate speed changes with r in accordance with the 
choice for a. In Fig. 1(a), the range for v grows steadily 
for a = 1, while the opposite is true for a = —1, as shown 
in Fig. 1 (b) . This outcome suggests that the eccentricity 
of the particle's orbit increases or decreases accordingly. 
When So/ [mor) = 10~^, it is clear from Figs. 1(c) and 
1(d) that the particle's motion becomes unstable due to 
the 0(^2) expression, and the co-ordinate speed rapidly 
approaches f = 1. For Fig. 1(c) with a = 1 this occurs at 
around t = 1500M, while Fig. 1(d) with a = —1 denotes 
this outcome at around r = 2000M. 

Since roughly the same behaviour occurs for the Kerr 
background in Paper I, this provides further evidc^ncc 
that a sufficiently large choice for Sq / {mo r) can trig- 
ger the transition from stable to imstablc orbital motion. 
This motivates a similar examination of the M0ller radius 
P(t) = (s/m)(r), as given by (32), and expressed by Fig- 
ure 2 for the same set of initial conditions, to third-order 
in e. When compared with Figure 1, and in general agree- 
ment with the corresponding plots in Paper I, there is 
further confirmation that the spin-curvature interaction 
due to p{t) induces the respective kinematic outcome 
for f (r). As well, the expression to third-order in e im- 
plies that the mass shift contribution fn"^ and the second- 
order spin shift term Sj listed in (32) results in a down- 
ward shift in the range of oscillation for Figure 2. Given 
so/{'mo r) = 10~2, comparison of Fig. 2(a) for a = 1, with 
Fig. 2(b) for a = —1 shows relatively little difference be- 
tween them. Not surprisingly, the situation changes con- 
siderably when So/ {mo r) = 10~^, as shown in Figs. 2(c) 
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and 2(d). It is very interesting to note that, when com- 
pared with Figs. 1(c) and 1(d), the sudden increase in v 
for both cases corresponds precisely with the condition 
that p{t) < for r > 1500M and r > 2000M, respec- 
tively. This provides more evidence suggesting that the 
M0ller radius must be strictly positive-valued to corre- 
spond with stable orbital motion, roughly agreeing with 
a similar set of conditions noted in Paper I for the Kerr 
background. 

Following the treatment given in Paper 1, the time- 
averaged value for the M0ller radius (p) = (s/m) as a 
function of the initial spin orientation angles 6 and is 
considered, leading to three-dimensional plots given by 
Figures 3 and 4 for so/{mor) = 10~^. In similar fashion 
to that shown in Paper I, these set of plots identify an 
even function symmetry according to = tt. As denoted 
by Figs. 3(a) and 4(a), there exists a non-trivial peak 
and valley structure in (p) that agrees in form with the 
corresponding set in Paper 1, while Figs. 3(b) and 4(b) 
indicate a loss of structure in < ^ < tt, with two peaks 
that remain. However, unlike the plots given in Paper I, 
this set of plots also numerically agree with each other, 
irrespective of the choice for a, which suggests that the 
relevant terms are either small compared to terms not 
coupled to a, or they integrate to zero entirely. 

As in Paper I, examination of the linear momentum 
components P^(r), based on (24) and (28), is in order. 
This is given by Figures 5-7, which show the radial, po- 
lar, and azimuthal components of the linear momentum 
in the Vaidya background, while the ratio P^{t)/P°{t) 
is displayed in Figure 8. Concerning the O(e^) expres- 
sion for so/(mor) = 10^^ and a = 1, Fig. 5(a) exhibits 
a slight contraction in the amplitude before outwardly 
expanding, with a similar behaviour shown in Fig. 5(b) 
for a = —1. Unlike the corresponding set of plots shown 
in Paper 1, the difference between these two plots is vir- 
tually negligible. When considering so/(mor) = 10"-^ 
for a = 1, the outward growth of the 0(£^) expres- 
sion for Fig. 5(c) starts to become dominant just before 
T = 1500M, in accordance with the growth of v{t) in 
Fig. 1(a). This is also true for Fig. 5(d) when a = — 1 at 
around r = 2000M. 

Figure 6 shows the polar component of the linear 
momentum, which indicates that while the expression 
to first-order in s remains around zero on average, 
the expression to second-order is positive-valued. For 
So/{mor) — 10^^, both Figs. 6(a) and 6(b) for a = 1 
and a = —1, respectively, display a slight positive- valued 
magnitude due to the 0(£^) expression for P^{t). This 
behaves similarly to the result obtained in Paper I for the 
Kerr background, also indicating that the spinning parti- 
cle will no longer remain on the orbital plane after a suffi- 
ciently long time. Furthermore, in terms of the choice for 
a there is essentially no difference in magnitude between 
these two plots for so/{mor) = 10~^. When considering 
so/{mor) = 10~^, as shown in Figs. 6(c) and 6(d), the 
expression to second-order in e is more obviously non- 
zero compared to the first-order contribution, indicative 



of the orbital instabilities suggested by Figs. 1(c) and 
1(d). 

The azimuthal component of the linear momentum is 
given by Figure 7, where Figs. 7(a) and 7(b) refer to 
a = 1 and a = —1, respectively, for so/(mor) ~ 10^^. 
In both cases, the O(e^) expression has very little im- 
pact on the overall plots. Furthermore, it is evident that 
the last term in (58d), which is linearly time-dependent, 
is most likely responsible for the plots' slopes, due to an 
overall minus sign for this term. When So/(too f) = 10~^, 
Figs. 7(c) and 7(d) indicate a strongly increasing ampli- 
tude for the expression to second-order in e, in similar 
fashion to the corresponding plots shown in Paper 1. 

The ratio P^(t)/P°(t) is given by Figure 8 for 
Sf)/{mf)r) = 10~^, where Fig. 8(a) describes a = 1 and 
Fig. 8(b) refers to a = — 1. Focussing on the expression to 
first-order in £, it is interesting to note that the amplitude 
shows a steady increase (a = 1) or decrease (a = — 1), 
unlike the strictly constant amplitude for the correspond- 
ing plots fomid in Paper I for the Kerr background. This 
is due to the last term of (59), which is linearly time- 
dependent with an overall positive sign. Similarly to the 
plots of Paper I, when adding the second-order contribu- 
tion, the amplitude for the ratio slightly contracts before 
steadily growing in magnitude. 

As a final example, it is useful to examine one of the 
components of the spin tensor to illustrate its properties 
due to spin-gravity interaction. In keeping with Paper 
I for the sake of comparison, the 5'°^(r) is chosen for 
study. This is given by Figure 9, expressed to third-order 
in £. Focussing on Figs. 9(a) and 9(b), corresponding to 
a = 1 and a = —1 for so/{mor) = 10^^, it is clear 
that the direction of radiation flow into or away from the 
black hole impacts upon the amplitude of S'"^(r) about 
zero, and that the expression to third-order in s deviates 
slightly from that due to the second-order contribution 
alone. This is in contrast to the corresponding set of plots 
in Paper I, which show no significant difference between 
the expressions to second- and third-order in e. When 
considering so/{mor) = 10~^, as shown in Figs. 9(c) 
and 9(d) for a = 1 and a = —1, respectively, the large 
outward growth of the amplitude due to the expressions 
to second- and third-order in £ is a further response to the 
orbital instability experienced by the spinning particle, as 
reflected in Figure 1. 



VI. CONCLUSION 

This paper is an application of the generalized CMP 
approximation approach to the Mathisson-Papapetrou- 
Dixon equations of motion of a spinning point particle 
in orbit around a spherical black hole in the presence of 
radially inflowing and outflowing radiation, as described 
by the Vaidya metric. When compared to a similar anal- 
ysis performed for orbital motion around a Kerr black 
hole, as described in Paper I [11], all relevant computa- 
tions, including the "radiative corrections" to the parti- 
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clc's squared mass and spin magnitudes, have nontriv- 
ial properties due to the exphcit time-dependence of the 
space-time background. It is somewhat ironic that, while 
the Vaidya metric is arguably much simpler in form than 
the Kerr metric, its time-dependence leads to much more 
complicated mathematical structure in the generalized 
CMP approximation than displayed in the previous ap- 
plication. As with Paper I, some numerical analysis is 
performed to illustrate conditions for the emergence of 
instabilities in the particle's orbit, with the suggestion 
that the M0ller radius needs to remain positive- valued in 
order to avoid the transition away from stable motion. 

The next step in this exploration is to obtain the per- 
turbed orbit from the results determined with the gener- 
alized CMP approximation, while also incorporating the 
effects of gravitational radiation within the process. As 
noted in Paper I, this generalization introduces concep- 
tual and technical challenges that are still not clearly 
understood at present. This consideration will be de- 
ferred to a future publication once these challenges are 
overcome. Another possibility is to explore a many-body 
interaction with spin incorporated through the general- 
ized CMP approximation. To do this requires under- 
standing the expected tidal and spin-spin interactions to 
be found when dealing with such a problem, which have 
a separate set of conceptual and technical challenges to 
consider. Nonetheless, both the work presented here and 
in Paper I illustrate the potential that comes from this 
line of research. 



which leads to 

A^S\l}-B, = 0, (A.l) 

where 

= -Pli^Sj°J. (A.2b) 

With (A.l), the S^yj^ components can be solved alge- 
braically in terms of the purely spatial components S^j^ , 
while the three spatial components are determined as so- 
lutions to the matrix differential equation 

^ = ^H-2-(o)r%.4^i = 0. (A.3) 
In explicit component form, (A.3) is equivalent to 
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dr 



+ 7;o^''SI>{t) = <5i2(t), (A.4a) 



^^ + lP'^S^}\r) = fe(r), (A.4b) 
^^ + ll''Sl^\r) = J3i(r), (A.4C) 



APPENDIX A: SOLUTION FOR THE 
FIRST-ORDER SPIN SHIFT 

To solve for sf directly from (26) requires use of the 
spin condition constraint equation (10) for j = 1 [11], 



where a*-^ , /3*-^ , and 7*-^ are antisymmetric spatial tensors, 
which may be r-dependent for a given choice of metric. 
For the Vaidya metric given by (38), with (45) chosen for 
AM and recalling (60), it is shown that 
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d(AM) 



a 
r 



dC 
d(AM) 



Okt a [1 — cos(2 0kt) 

NA^ ~ m 



r Op 



d? 



^ - 2 (r^ O^)] [1 - cos(2 Ok t)] | 
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„23 



and 



(A.5a) 
(A.5b) 

(A.5c) 
(A.5d) 

(A.5e) 
(A.5f) 



S23{t) 



1 mo r 



A N'^A \mor 
d(AM) 



a 
A 



d? 



(r^O^) 3 Q,- (Okt, 2^,,^) 



— (r^O^,) Q-(3OKr,20,0)-3g-(OKT,2^?,-^ 



- ^(rOK)Qr(^^KT,2 0,.^)(O 



kt)| 



3 mor^ 



2 



d(AM) 



d^ 



So 
mor 



(r" O2,)Q7(OKr,2 0,(/.) (Ok r) , 



1 mor f So 



\mor 
d(AM) 



]4 



dC 



(r^Ol) 30+(2OKr,2^,2(^) -3 sin(20Kr-2^) 



(r^ 02,) [q+ (4 Ok r , 2 ^ , 2 (^) - 2 cos(4 Ok r - 2 



+ 12Q+(2OKr,2(9,0) + 8 cos(20Kr) + 3Q+(0,2^,2<^) -6 cos(2^) 



+ ;^(rOK) g3+(20KT,2^,2(/.)-2sin(20Kr-2<^) (Okt) 



r 



(A.6a) 
(A.6b) 



(A.6c) 



which are substituted into (A. 4) to solve for S'^1?(t). 



APPENDIX B: ORTHONORMAL TETRAD 
FRAME 



thonormal tetrad frame A^^ corresponding to fundamen- 
tal static observers in the Vaidya background, subject to 



This Appendix outlines the derivation of the or- 

thonormal tetrad frame A^q, for orbital motion in the 
Vaidya space-time background, following the approach 
given elsewhere [8, 22], with the assumption that 
AM/Mq < 1. To proceed, recall from (38) the Vaidya 
metric g^i, in {t, r, 9, 0) co-ordinates and consider the or- 



Va0 = 5/x,.A^sA' 



/3 • 



(B.l) 
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The tetrad set for static observers is assumed to take the 
form 



(Zo, 0,0,0), 
(Zi,Z,0,0), 



= (0,0,-, 



A^o = 0,0,0, 



r sm( 



where satisfying (B.l) leads to 



Z = 
Zo = 
Zi = 



A' 



2 a AM 



A^Z 



(B.2) 
(B.3) 

(B.4) 
(B.5) 



(B.6) 
(B.7) 
(B.8) 



The main idea is to Lorentz boost A'*^ — » A'*^ with 
speed /3 to its location on the orbit, set A^q = A^q, and 
determine A'^j accordingly to accommodate the parallel 
transport condition DA'^a/dr = 0. Because M changes 
with proper time r along the null co-ordinate ^ according 
to (37), it follows that the frame will be boosted along 
at least the azimuthal and radial directions. This implies 
that the orbit is strictly no longer circular and introduces 
some complications in determining A''q. However, it is 
shown in Appendix C that, for ^ = 7r/2, 



(B.9) 



where 



K 



4 a AAf 



L 



+ r Au^ 



where 7 = l/yl — is the Lorentz factor. It follows 
that identification of (B.ll) with (B.9) leads to 

Z 



7 = 



/3 = 



K 



E ^ r. 2a AM , ; 

— +Au°- — — Au^ 

A^ A^ Z^ r 



(B.15) 



(Auy 
Z2 



r Am' 



r 



1/2 



a 



tan 



L , A«i 
-+rAu^] —— 
r Z 



, (B.16) 



(B.17) 



Finally, the spatial triad need to be rotated back by 
r^KT to reflect the parallel propagation of the tetrad 
along the orbit [22], such that 

A'^i = A**! cos(Okt) - A'^g sin(OKT), (B.18) 
A'^2, (B.19) 



A^, 



A'^i sin (r^K t) + A^'g cos (Ok t) . (B.20) 



This leads to the final expression for the orthonormal 
tetrad frame 



A''. 



^(| + A.o,A.S0,i(^ + .A.3)), 

(B.21) 



A^ = 



-7 ^ sin (Ok t) + ^ ^ (Ok t, a) 



ZF+{nKT,a),0,-Ff (Ok r, a) ) , (B.22) 



0,0, -,0 
r 



(B.23) 



7 P cos (Ok r) + ^ ^ (Ok t, a) 



ZF+iflK T, a) , , - J ^TC^^K r, a)^ , (B.24) 



1/2 



(B.IO) 



where 

F±(OKr,a) 



is a normahzation condition for A^'g, and Att°, Au^, and 
Au^ are contributions to the tetrad frame's overall four- 
velocity determined in Appendix C. It is straightforward 
to confirm that K ^ 1 and A''g 
0, as expected. 

Applying a Lorentz transformation to A^a and given 
(B.9), it is true that 



A^Ocsch) whenAM 



- (1 - 7) cos (Okt - a) 
±^(l + 7)cos(OKT + a), (B.25) 
F±(OKT,a) = ^(1-7) sin (Okt -a) 

±i(l+7)sin(0KT + a). (B.26) 

As a final consistency check, it is straightforward to ver- 
ify that A^3 (Ok r + f ) = A'^j (Ok r), and that (B.21)- 
(B.24) reduce to (40a)-(40d) in the limit as AM 0. 



A^g = 7 A'"o + /3(cosaA''3 +sinaA''j) 
A'^j = cosa A^j - sina A^3 , 
A'^S = A'^s, 



A'^g = 7 (cosaA'^3+sin(5A^j)+/3A^o 



(B.ll) 

(B.12) 
(B.13) 
(B.14) 



APPENDIX C: ZEROTH COMPONENT OF THE 
TETRAD FRAME 

This Appendix outlines the method to determine the 
zeroth component A'^g of the orthonormal tetrad frame in 
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the Vaidya space-time background. Suppose that X^q = 

K^q/K, where 



6 (Sch) 



AA^, 



AA'^o = (Au°,Aw\Au2,Au3) , 



and 



1/2 



(C.l) 

(C.2) 
(C.3) 



Prom (38), it is possible to identify the metric connection 



as 



(C.4) 



where AV^ocf) represents the contributions dependent on 
AM. Then it foUows from requiring DA'^g/dr = that 



d 
d^ 



- (AA%) 



P^„(AA"o) = 



(C.5) 



where 



= -(ArV) A"6(Sch)A^ 



(Sch) ^ (Sch) ' 



(C.6) 
(C.7) 



After specifying 6* = 7r/2 for the orbital plane, it is shown 
that d (AA^q) /dr = = o, which impHes that AA^q 
is a constant that can be set to zero. This leads to the 
column vector differential equation 



/ AAOq \ 









2E Mo 





\ 


/AA% 


\ 




AA^o 


+ 




2EMn 





2A^ L 
r 




AA^o 




= Q' 


\ AA^Q ) 




I 





2L 





/ 


Vaa^o 


/ 


Uv 



to solve, where 



for 



Q° = 



d (AM) 



/AMY 



2L2 am 



= 0. 

In matrix notation, (C.8) is represented as 
d 



dr 



(AAo) + 7'(AAa) = Q. 



(C.9) 

(C.IO) 
(C.ll) 

(C.12) 



The normal mode expression for AA5 is determined by 
finding an invertible matrix C that is constant in r, such 

that for AA^ = C'^ AAg and q' = C'^ Q, 



^(AA1)+P' (AA^) 



- / 



(CIS) 



where V = C ^VC is diagonalized. It follows that 
^00 ^ \ 



V' 



2i{rQ.K) 
\0 -2i(rOK)/ 



\ 1 



E{rn-K) ■ E{ra-K) \ 
—43— 43— \ 



, (C.15a) 



/ EL a E^jr^K) \ 
I -pr U 33 — \ 



_ . E(rnK) 1 ■ A^L 
'2 2 ' 2(rf2K) 



g(rOK) 



2 * 2(rnK) / 



(C.14) leading to 



(C.15b) 
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_d_ 

dr 



AV" 











2i{rflK) 
V 



I 



(C.16) 



where 



Q" = 

The solution to (C.16) is then 



AX 



rO 



Q'%r')dT'., 



(C.17) 

(C.18) 
(C.19) 

(C.20) 



AA'^Q = e-2^("«^) re'<"«"')Q'i(r')dT', (C.21) 
Jo 

(C.22) 



AX% = {AX%y , 
which from AAq = C AAg and (C.2) leads to 

A«o 



E {rftKY 



Au'^ = 2Re(AA'^Q) , 
At(2 = 0, 



Aw^ = AA'"fi + 



2L 



° r2 (rOK) 



lm(AA'i6) . 



(C.23) 
(C.24) 
(C.25) 

(C.26) 



the nonzero components of the Riemann curvature tensor 
^RjXi>ce$ Fermi frame are listed as follows: 



APPENDIX D: FERMI-FRAME RIEMANN 
TENSOR COMPONENTS 



Given that the nonzero Riemann tensor components 
for the Vaidya metric are 



Rowi = --^{Mo + AM). 
i^0202 = — (Mo + AM) H 



R 



d(AM) 



a 
l2 



-R0212 — —tt; 
-R0313 



-2 [Mq + AM) 



AM 



d (AM) 



-R0202 sin 9 , 
i?02i2 sin^ 9 , 
1 



1212 — 



^4 



1 / 2AM\ 
r V r 



j (Mo + AM) + 



(D.l) 
(D.2) 

(D.3) 

(D.4) 
(D.5) 
d(AM) 



-R1313 = -R1212 sin 9 , 

-R2323 = 2r(Mo + AM)sin2 6', 



(D.6) 
(D.7) 
(D.8) 



20 



■0101 



(Au^)^ iioiOl + + ^0303 



sin^ (Ok t) 



+ 



2/37 

a:2 



£; ^ n 2 a AM ^ , , ^ , „ 



L ^ .V ( r. '^Oi AM ^ 

+ AuM i?0313 + Ri 



^2^2 J. 



0303 



F+(OKT,a) 



Au^ i?0313 +(^+Au°] i?0303 



(Qk r, a) \ sin (Ok r) 



I f E ^0 2a AM ^ i\ „ 



^ + Am'' 



4 / 1 A7\ 

^1313 + -^r^ i?0313 + — 



1 AM \ AM 



^0303 



[F+(OKT,a)]' 



2Z /L ^3 
^ + Au^ 



\ , 2a f E ^ A AM „ 

AU^ iil313 + -To-T^ + Au° i?0303 



y42 ^2 V ^2 



S ^ n 2 a AM ^ , , „ 

]42+^« +]42:^ — ^" 1^0313 



j;-(OKT,a) F+(OKr,a) 



^{(A.^)^i?i3i3+(^+A.o 



— + Am° ) i?0303 + 2 Au^ i?0313 



[F-(OKr,a)]' , (D.9) 
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0103 



■0101 



4 + Au^ 



0303 



sin^ (Ok t) cos (Ok t) 



iC2 \ 



E ^ n 2 a AM ^ , , ^ , 



^ . 3\ /^r, 2a AM „ 
^ + A«3 i?03i3 + B, 



■0303 



Au^ i?0313 + + Au°^ i?0303 



F-(2nKT,a) 



A2 



+ Au^ 



Z2 



1 AM 



AM 



Au^ 



0101 



L 



+ Am' 



4 / 1 

^1313 + ^2 ^2 I "^0313 + ^2 ^2 



1 AM \ AM 
- -^0303 I 



F+inKT,a) F+(r!KT,a) 



^ + Au^ 



^2 J. I 



\ , 2a f E ^ n\ AM „ 

A«l i?1313 + f + Au° j — i?0303 



, .0 2a AM . n „ 

+ ( + A^O + ^ Au' ] i?0313 



[F- (l^K T, a) (l^K r, a) - i^" (f^K r, a) (I^k t, a) 



- { M ' ^1313 + + Aii' 



+ Am° ) i?0303 + 2 Au^ i?0313 



xf-(nKr,a)j;-(nKr,a) , 



(D.IO) 
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I Au^ iJoioi F+ (0, a) + ^(^^+ Au^^ i?0303 Ff (0, a)| sin (Qk r) 



pjZ f E ^0 2a \ 



[F+ (Ok r, a)] ^ cos (Ok t) + F+ (Ok r, a) sin (Ok r) 



/37 



]42 



+ Aw° 1 Eosos + Awl i?0313 



f-(0,a) F-(OKr,a) 



+ 4^ ( 4 + Au3 ) ( i?03i3 + ^ i?0303 ) [f^- (2 Ok r, a) (Ok r, a) 



y42 ^2 



2 F- (Ok r, a) F+ (Ok t, q) sin (Ok t) 



a: 



z r /L ^ o\ r 4 / 1 ATI 

[^[■^ + ^'' ) [-^1313 + ^0313 + — 



1 AM ^ \ AM 



^0303 



F+ (Ok r, a) 



\ 1 „ 2a f E ^ o\ AM 

A«l i?1313 + (^^ + Ah" j — 



f E ^ 2a AM ^ A „ 

i?0303 + (^]42 + + — ^« j ^0313 



X F- (Ok r, a)} [F" (Ok t, a) F+ (Ok t, a) + F" (Ok r, 5) F+ (Ok t, a)] , 



(D.U) 



— + Aw° ) i?0202 + 2 Au^ i?0212 



^+Aw3^ i?2323| (D.12) 



0212 



/37 



E 

-^+Au°] Rq202 + Au^ i?0212 



sin (Ok r) + ( ^ + AuM E2323 F" (Ok t, a) 



Au^ -R1212 



2a f E 



, n \ AM f E , n 2 a AM , , 

Aw" ] i?0202 + ( — + Au" + 



Au^ iJ, 



1-0212 



X F+(OKr,a) , 



(D.13) 



0223 



01 

KZr^ 



E 



+ Am" I i?0202 + Av} FiQ2\2 



1 L 



cos (Ok r) + ( ^ + Au^ ) i?2323 F^ (Ok t, a) 



. 1 „ 2 a / F , AM „ / F , n 2 a AM , , , „ 

AW^ i?1212 + [-^ + A^^" j — iio202 +(^^42+ A«° + — A^xl ) i?0212 



xF+(OKr,a) , 



(D.14) 
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0303 



(Au^)^ i?0101 + + ^0303 



cos^ (Ok t) 



E , n 2 a AM ^ , , ^ , „ 



L A 3^ 2a AM 

— + Au I ( -Ro313 + ^2 ^2 -^0303 



F+(OKT,a) 



1 / L 



Zr Kr"^ 



+ Am' 



Am^ Eo313 + (^-J2+ ^0303 



^7 (2 Ok r, a) ^ cos (Ok t) 



Z2 / S . n 2 a AM , , , 



0101 



Au- 



4 / 1 

-^1313 + ^2 I " -^0313 + ^2 _^2 



1 AM ^ \ AM 



-Ro303 



[i^+(OKr,a)] 



2Z /L ^3 
^ + Am^ 



J^2 J. (^ J. 



. 1 „ '2a / E , n\ AM „ 

Am^ iil313 + f + Am° j ^ i?0303 



, ^ n 2 a AM ^ , , „ 

+ ( + AU° + — Au' ] i?0313 



j;-(OKT,a) i?+(OKT,a) 



+ 



^|(Am1)'J?1313+(^+Au°) 



+ Au° I i?0303 + 2 Am^ i?0313 



[F-(OKr,a)]' , (D.15) 
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■0313 



KZ 



I Aui iJoioi F+ (0, ") + + ^"^) ^0303 -Fs" (0> a)} cos {SIk t) 



l3jZ ( E 



K 



2a 



— + Au° - 1 iioioi 



[F+ (Qk t, a)] ' sin (Ok t) + F+ (Ok r, a) cos (Ok t) 



+ 



/37 



+ i?0303 + A^i i?, 



10313 



J^-(0,a) i7(0KT,a) 



/37 



+ I — + Au^ ) i?0313 



^ " it!o303 ) [J^c+ (2 Ok r, a) i^" (Ok t, a) 



^2 



+ 2 F- (Ok r, a) F+ (Ok r, a) cos (Ok t) 



^ ^Zr4 + A«^^ 



4 / 1 

^1313 + ^2 ( " -^0313 + ^^5-^ 



1 AM ^ \ AM 



-^0303 



F+(OKr,a) 



1 

+ - 

r 



■ i„ 2a [ E ^ o\ AM „ /F ^ q 2q AM ^ 

A«l iil313 + (^^ + A«OJ iio303 + (^]42 + + ^42^ 1 ^0313 



F- (Ok t, a)] [F' (Ok r, a) F+ (Ok r, a) + F' (Ok t, a) F+ (Ok t, a)] , 



(D.16) 



R 



1212 



/?'7' „ . 2,0 ^ 2/37 2a AM 

" I -K0212 



Z2l 



-R0202 sin (Ok t) 



+ 



A2Z2 r 



^0202 I F+ (Ok t, a) sin (Ok t) 



+ ■ 



4 / „ 1 AM „ \ AM 

-R12I2 + r,0 I -R02I2 + .» „o -R02O2 



A2 Z2 



^2^2 r 



[F+(OKT,a)] 



+ -^^^2323 [F+(OKr,a)]' , 



(D.17) 



R 



1223 — 2^2 r2 



^ ^ -R0202 sin (Ok t) cos (Ok t) - — ^ f -R0212 + .o"; i?0202 ) -Fc"*" (2 ^^k t, a) 



Z2 



Z2 
^2 



-R1212 + 



A2 Z2 



1 AM \ AM 

OL -R0212 + ^2 ^2 ^ -R0202 



F+(Okt,q) F+(OKr,Q) 



+ 4: ^2323 -Fs" (f^K r, 5) F- (Ok r, a) 



(D.18) 
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^1313 = 7' { ^ioioi [ ( [F+ {Qk t, a)] ' cos^ (Ok r) + [F+ (Ok r, a)] ' sin^ (Ok r)) 



+ 2 F+ (Ok t, a) F+ (Ok t, a) sin (Ok r)] 



+ ;^272 -^303 



( [F- (Ok t, a)] ' sin^ (Ok r) + [F" (Ok r, a)] ' cos^ (Ok t)) 



2 F- (Ok t, a) F" (Ok r, a) sin (Ok r)] } 



2 3', . „ 
H 5— -K0313 + 



2 a AM 



A^Z^ r 



i?0303 i^s" (0, a) [F- (Ok t, a) F+ (Ok r, a) + F" (Ok t, a) F+ (Ok t, a)] 



Z2 
^2 



-R1313 + ^2 ^2 ( '^-^0313 + 



1 AM „ \ AM 



^2^2 ^ -R0303 ^ 



-sn2 



[F- (Ok r, a) F+ (Ok t, a) + F" (Ok r, a) F+ (Ok r, a)] 



(D.19) 



R 



2323 — ^2^2 



^ ^ -R0202 cos^ (Ok t) + ^i?02i2 + ^2^2 ^0202^ F+ (Ok t, a) cos (Ok t) 



+ — 



1212 



a Ri 



1 AM „ \ AM 



•0212 



712^2 r 



Rq 



202 



[F+(OKr,a)] 



+ ;^ii2323 [F-(OKr,a)]' . 



(D.20) 
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FIG. 1: Co-ordinate speed u(r) of the spinning particle for circular motion around a black hole described by the Vaidya metric, 
for r = 6M and 9 = cj> = 7r/4. Fig. 1(a) shows a gradual increase in amplitude for so/{mor) = 10"'^ with no significant 
difference due to second-order contributions in e. For a = —1, Fig. 1(b) shows a corresponding gradual decrease in amplitude. 
In contrast, Figs. 1(c) and 1(d) indicate the existence of an instability in the orbit when so/{mor) — 10^^, due to the rise of 
v{t) from the second-order contribution in e. 
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FIG. 2: M0ller radius p(t) = {s/m){T) in the Vaidya background for r = 6M and 6 = (j> = 7r/4, in units of so/mo. Figs. 2(a) 
and 2(b) show that while the higher-order contributions in e lead to a slowly increasing amplitude in p, Figs. 2(c) and 2(d) 
indicate a much larger amplitude increase as so/(mor) — 10^^, where the second- and third-order contributions in e become 
distinctive. 




FIG. 3: Three-dimensional plot of 
the time-averaged M0ller radius (p) — 
(s/m) as a function of 9 and </) for 
r = 6M and a = 1. Fig. 3(a) shows 
a complicated peak and valley struc- 
ture to (p) that simplifies somewhat 
in Fig. 3(b). 



(a) a = l,0{e^) (b) a = 1 , 0{e^) 




FIG. 4: Time-averaged M0ller radius 
as a function of 9 and (j) for r = 6A/ 
and Q = —1. The plots are essentially 
indistinguishable when compared to 
Figure 3, as a = 1 goes to a = — 1. 
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FIG. 5: Radial component -P^(t) of the linear momentum in the Vaidya background for r = 6M and 6 = <f) — 7r/4. When 
adding the second-order contribution in e for so /{mo r) = 10~^, Figs. 5(a) and 5(b) indicate a shght increase in the amphtude, 
while Figs. 5(c) and 5(d) show a much stronger amplitude increase as so/{mor) — 10~^. 
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FIG. 6: Polar component P^{t) of the linear momentum for r = QM and 6 — (j) = tt/A: in the Vaidya background. All three plots 
show that the second-order contribution in e results in a net non-zero magnitude for the polar component. While Figs. 6(a) 
and 6(b) yield a modest non-zero effect for so/(mor) — 10~^, Figs. 6(c) and 6(d) show a significantly more pronounced effect 
for so/{mor) — 10~^ . 
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FIG. 7: Azimuthal component P'^{t) of the linear momentum for r — QM and 6 — (f> = tv/A in the Vaidya background. 
For Figs. 7(a) and 7(b), the change in magnitude is almost totally dominated by the first-order contribution in e. Figs. 7(c) 
and 7(d), in contrast, indicate that the second-order contribution in e dominates with a strongly unbounded increase in the 
amplitude as so /{mo r) — 10~^. 




FIG. 8: Ratio of P^(r) to P°(t) in the Vaidya background for r = 6M and 9 = (f> = n/A. Figs. 8(a) and 8(b) indicate that the 
second-order expression in e leads to a wider variation in the ratio than found in the first-order expression only, irrespective of 
the given choice for a. 
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FIG. 9: The S'^'^ir) component of the spin tensor for r = 6Af and 6 



7r/4 in the Vaidya background. Fig. 9(a) shows a 



modest decrease in amplitude for so /{mo r) = 10 and q = 1, with a corresponding modest increase in amplitude in Fig. 9(b) 
for a — —\. In contrast, both Figs. 9(c) and 9(d) show a significantly more pronounced amplitude increase for so/{mo r) = 10"^ . 



